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FOREWORD 

The  work  described  in  this  report  was  carried  out 
as  a  part  of  a  program  to  investigate  the  effects  of 
mass  transfer  on  the  stability  of  re-entry  type  vehi¬ 
cles.  In  addition  to  the  work  reported  here,  theoret¬ 
ical  investigations  were  made  of  the  aerodynamic  proc¬ 
esses  involved  and  these  are  described  in  a  companion 
report.  The  overall  program  was  supported  by  the  Ad¬ 
vanced  Research  Projects  Agency,  Ballistic  Missile 
Defense  Office,  Penetration  Aids  Branch,  and  was  tech¬ 
nically  administered  by  the  Fluid  Dynamics  Branch  of 
the  Office  of  Naval  Research.  The  development  of  an 
apparatus  as  complex  and  sophisticated  as  that  em¬ 
ployed  in  these  studies  requires  the  expert  skills 
and  devotion  of  many  people.  The  authors  wish  to  ac¬ 
knowledge  the  contribution  of  their  co-workers  in  gen¬ 
eral  and  in  particular,  Messrs.  DeCoursin  and  Lindahl 
for  contributions  to  the  development  of  the  concepts 
of  test  apparatus,  and  Messrs.  Estel  and  Matsuura  for 
their  contribution?  to  the  solution  of  the  many  de¬ 
tailed  problems  which  arose  in  converting  these  con¬ 
cepts  i ntc  a  working  apparatus. 
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ABSTRACT 

This  report  describes  an  experimental  study  of  the 
effects  of  mass  transfer  on  the  stability  of  conical 
bodies.  A  unique  apparatus  was  created  which  permits 
the  support  of  a  conical  model  in  a  hypersonic  wind  tun¬ 
nel  with  small  pitch  damping  due  to  the  apparatus,  and 
at  the  same  time  provides  a  means  for  blowing  through 
four  separate  (aft,  forward,  top,  bottom  combinations) 
model  surface  areas.  Independent  servo  controls  were 
provided  for  the  flow  through  each  surface  section  such 
that  the  mass  transfer  through  the  model  wall  could  be 
varied  both  spatially  and  temporally.  This  apparatus 
was  installed  in  the  FluiDyne  Hypersonic  Wind  Tunnel 
and  tests  were  run  at  a  nominal  Mach  number  of  11.  Runs 
were  made  without  mass  transfer,  with  mass  transfer  pro¬ 
duced  by  a  subliming  material,  and  with  mass  transfer 
produced  by  blowing  through  a  porous  surface.  These 
tests  provided  additional  confirmation  that  mass  trans¬ 
fer  from  a  surface  will,  in  itself,  significantly  in¬ 
fluence  the  pitch  damping  characteristics  of  a  conical 
shape,  whether  its  origin  is  in  a  heat  transfer  related 
phenomenon  (such  as  ablation  or  sublimation)  or  an  in¬ 
dependently  controlled  flow  process  such  as  is  repre¬ 
sented  by  blowing  through  the  conical  surface.  In 
addition,  these  tests  provided  an  opportunity  for  shake- 
down  of  the  blowing  flow  control  apparatus;  however,  it 
was  not  possible  to  carry  the  blowing  study  into  the 
parametric  investigation  initially  planned  within  the 
present  program.  Thus,  the  interrelation  of  such  para¬ 
meters  as  mass  flow,  volume  flow,  spatial  distribution, 
temporal  distribution,  etc.,  remain  to  be  investigated 
in  detail. 
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FiuiDyne  Engineering  Corporation 

1.0  INTRODUCTION 

This  report  describes  the  experimental  apparatus 
which  was  developed  to  study  the  influence  of  mass 
transfer  from  a  surface  on  the  dynamic  stability  of  a 
conical  vehicle  during  hypersonic  flight,  and  the  re¬ 
sults  of  a  series  of  wind  tunnel  tests  employing  this 
apparatus . 

Exploratory  wind  tunnel  tests  which  were  conducted 
at  FluiDyne  Engineering  Corporation  in  1963  demonstrated 
that  mass  transfer  from  the  surface  of  a  conical  re-entry 
body  can  have  a  significant  effect  on  the  pitch  damping 
coefficients  of  such  a  body  (Reference  1).  The  magnitude 
of  this  effect  was  found  to  be  sufficient  to  reverse  the 
natural  stability  characteristics  of  the  body,  provided 
the  distribution  of  the  mass  transfer  was  of  the  necessary 
form.  Very  generally,  mass  transfer  from  areas  forward 
of  the  center  of  oscillation  appeared  to  stabilize  a  booy 
whereas  mass  transfer  from  areas  aft  of  the  center  of 
oscillation  tend  to  destabilize  it.  The  importance  of 
these  effects  as  regards  vehicle  design  and  the  prediction 
of  performance  character i st i cs  was  immediately  apparent. 
Thus  the  existence  of  a  significant  aerodynamic  phenomenon 
had  been  demonstrated;  the  immediate  problem  then  became 
one  of  analyzing  this  phenomenon  in  sufficient  detail, 
by  both  theoretical  and  experimental  means,  to  identify 
the  governing  parameters  involved  and  develop  the  equations 
relating  these  parameters  to  the  observed  results. 

In  the  case  of  actual  flight  vehicles  the  basic 
mechanism  by  which  mass  is  transferred  from  the  vehicle 
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surface  is  the  ablation  of  material  placed  on  the  surface 
to  protect  the  interior  from  overheating.  From  both  a 
theoretical  and  experimental  point  of  view,  i.he  inclusion 
of  parameters  which  define  the  ablation  process  itself 
into  the  study  of  the  effects  of  the  product  thereof  on 
vehicle  aerodynamics  leads  to  an  extremely  complex  pro¬ 
blem.  Hence,  while  it  was  recognized  that  the  ablation 
process  and  the  aerodynamic  flow  process  are  interrelated, 
it  seemed  desirable  to  attack  the  experimental  study  cf 
the  aerodynamic  process  as  a  mass  transfer  problem  without 
attempting  true  simulation  of  the  ablation  process.  Although 
the  aerodynamic  process  itself  is  very  complex  (the  effects 
of  temporal  and  spatial  variations  in  mass  addition  and 
the  effects  of  the  properties  of  the  fluid  being  added 
must  be  expected  to  be  significant  and  specifically 
accounted  for)  it  should  be  possible  to  develop  a  basis 
for  predicting  the  effect  on  the  overall  flow  process, 
and  thus  the  vehicle  dynamics,  of  a  given  mass  addition 
flow.  With  this  in  hand  the  characteristics  of  the  ablation 
process  which  interact  significantly  with  the  1  low  process 
as  regards  vehicle  stability  can  be  identified,  and  the 
effect  of  a  given  ablation  process  predated. 

The  current  program  evolved  from  these  general  con¬ 
siderations,  with  early  experimental  work  being  performed 
on  a  side  supported  model  for  which  the  bearing  was  located 
outside  of  the  wind  +unnel  flow.  While  the  geometrical 
arrangement  in  this  early  apparatus  limited  the  degree  of 
simulation  which  could  be  achieved,  tests  were  run  with 
blowing  through  a  porous  surface  of  the  model  to  simulate 
ablation  and  the  results  obtained  (Reference  2)  conclusively 
established  that  the  basic  phenomenon  previously  observed 
as  a  characteristic  of  ablating  vehicle0  could  be  produced 
by  a  relatively  simple  mass  addition  flow. 
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Concurrent  with  this  early  experimental  work  theoretical 
studies  were-  carried  out  (and  reported  in  Reference  3)  which 
also  confirmed  that  the  observed  effects  could  be  produced 
by  a  simple  mass  addition  flow.  Both  of  these  early  re¬ 
sults,  while  genera' ly  significant  in  establishing  that 
a  complex  interaction  between  an  ablation  process  and  the 
flow  process  is  not  a  requirement  for  producing  an  effect 
on  vehicle  stability,  left  the  development  and  con.irmation 
of  a  means  of  predicting  performance  in  any  gi-'en  case  to 
be  achieved. 

In  order  to  provide  confirmation  of  a  theoretically 
developed  technique  for  predicting  performance  and/or  to 
provide  the  empirical  results  and  coefficients  often 
necessary  to  establish  a  base  from  which  a  basically  sound 
theory  can  be  used  to  predict  quantitative  results,  it  is 
necessary  to  carry  out  an  orderly  set  of  experiments  in 
which  the  various  major  parameters  can  be  varied  independent¬ 
ly  over  an  adequate  range.  In  the  present  case  this  require¬ 
ment  lea  to  the  design  an.J  fabrication  of  a  test  apparatus 
specifically  for  the  purpose  of  studying  the  effects  of 
varying  the  principal  mass  flow  parameters  on  the  stability 
of  a  conical  vehicle.  This  apparatus,  which  is  described  in 
detail  in  the  following  report,  consists  of  a  conical 
model  supported  byan  internal  air  bearing  w'ich  in  turn 
is  supported  from  a  conventional  aft  sting.  The  model 
is  compartmerited  so  that  blowing  forward  aft,  windward, 
and  leewa,'G  can  be  independently  varied,  and  the  flow 
through  each  compartment  is  controlled  by  a  servo  system 
such  that  the  flow  can  be  varied  with  time  in  terms  of 
the  form  of  flow  versus  time  curve  and  the  displacement 
of  this  curve  with  respect  t.o  the  model  motion.  Pro¬ 
vision  was  made  in  the  development  of  this  system  to  permit 


3 


FluiDyne  Engineering  Corporation 


its  use  in  testing  with  subliming  materials,  su*i  as 
par  ad i ch lorobenzene  and  ammonium  chloride,  and  with 
blowing  gases  other  than  air.  Because  of  the  complex¬ 
ity  of  the  sys.em.  a  great  deal  of  difficulty  was  en¬ 
countered  in  p.acing  it  in  operational  status;  thus  the 
air  bearing  with  its  integral  blowing  channels,  the 
frictionless  angle  of  attack  transducer,  the  minia¬ 
turized  flow  control  valves  and  their  e 1 ectrohydrau 1 i c 
servo  control  systems,  all  represented  the  development 
of  new  operational  systems  involving  sophisticated  en¬ 
gineering  concepts.  The  development  of  this  equipment 
culminated  in  its  use  in  carrying  out  additional  wind 
tunnel  studies  f  the  effect  of  mass  transfer  on  vehicl 
stability.  Tests  were  made  using  the  system  without 
mass  transfc  ,  with  mass  transfer  resulting  from  a  sub¬ 
liming  surface,  and  with  mass  transfer  produced  by  blow 
i ng  through  a  porous  surface.  The  results  of  these 
tests,  which  were  made  in  the  FluiDyne  Hypersonic  Wind 
Tunnel  at  a  nominal  Mach  number  cf  11,  are  described  in 
this  report.  In  addition  to  the  data  presented  here, 
schlieren  movies  of  several  of  the  test  runs  were  made. 
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2.0  FACILITY  DESCRIPTION 

The  tests  were  performed  in  the  FluiDyne  Hypersonic 
Wind  Tunnel  at  a  nominal  Mach  number  of  11.  This  is  a 
free  jet  tunnel  with  a  nozzle  exit  diameter  of  20  inches. 

Air  is  stored  at  5000  psi  and  is  manually  throttled  to 
the  desired  stagnation  pressure.  The  air  is  then  heated 
by  passing  it  through  a  zirconia  cored  brick  bed  storage 
heater  which  has  been  preheated  by  an  oxygen- propane  burner. 
A  sketch  of  this  tunnel  is  shown  in  Figure  1.  This  facil¬ 
ity  is  more  fuHy  described  in  a  booklet  entitled,  "Hyper¬ 
sonic  Aerodynamic  Studies,"  which  is  available  from  FluiDyne 
Engineering  Corporation. 

Stagnat'on  pressure  is  measured  with  a  test  quality 
bourdon  tube  type  gage  and  stagnation  temperature  with  a 
Pt/Pt-10  Rh  thermocouple  located  in  the  stilling  chamber. 

Run  times  are  on  the  order  of  60  seconds.  The  model  sup¬ 
port  system  provides  injection  and  retraction  from  the 
stream,  axial  translation,  and  pitching  from  -45°  to  +45° 
angle  of  attack. 

During  the  first  series  of  tests  (aolating)  in  June 
1966  and  during  the  second  series  of  tests  (blowing)  in 
January  1967,  nozz’e  calibrations  were  performed.  Results 
of  the  cal  ibrations  are  as  follows:  the  early  calibration 
gave  a  variation  of  +  0.8%  from  tne  mean  Mach  number  over 
an  eight  inch  test  core;  the  latter  calibration,  after 
nozzle  reassembly,  gave  a  variation  of  +  0.4%.  The  nozzle 
flow  characteristics  are  shown  in  Figure  2  in  terms  of 
test  :,‘ction  Mach  numbers  as  a  function  of  stagnation 
pressures  and  stagnation  temperatures. 
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3.0  TEST  MODEL  AND  INSTRUMENTATION  DESCRIPTION 

The  dynamic  stability  coefficients  for  both  ablat¬ 
ing  and  blowing  models  of  the  same  geometry  were  deter¬ 
mined  from  the  free  oscillations  on  an  air  bearing  support. 

In  the  case  of  the  blowing  model,  additional  ' nstrumentat i on 
controlled  the  blowing  so  as  to  simulate  ablation.  The 
subsystems  are  the  air  bearing,  the  model  position  trans¬ 
ducer,  the  blowing  flow  modulating  needle  valves,  the 
Pressure  transducers,  the  mass  transfer  control,  and  the 
data  acquisition  system. 

3.  1  Dynamic  Stability  System 

An  assembly  drawing  of  this  system  with  two  sectional 
views  is  shown  in  Figure  3A;  details  are  shown  in  Figures 
3B  and  3C.  A  conical  model  mounted  on  an  air  bearing  is 
free  to  oscillate  in  the  vertical  plane  about  its  center 
of  gravity.  In  the  case  of  a  blowing  model  as  many  as 
four  porous  skin  sections  may  be  employed. 

The  principal  element  of  this  system  is  the  air  bear¬ 
ing  itself.  The  design  of  an  air  bearing  which  would  oper¬ 
ate  satisfactorily  and  could  still  be  entirely  buried  within 
a  model  was  in  itself  a  major  problem  which  was  significantly 
complicated  by  the  need  to  "cross  the  bearing"  with  flow 
passages  to  tne  blowing  compartments  when  blowing  models 
were  employed.  Every  effort  was  made  to  assure  that  the 
bearing  oscillated  with  a  minimum  of  viscous  restraint  and 
with  no  mechanical  restraint.  To  accomplish  this  the  flow 
passages  were  machined  into  the  bearing  and  housing  elements 
using  what  amounted  to  labyrinth  seals  at  the  intersections 
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to  limit  leakage  to  insignificant  amounts,  and  no  electri¬ 
cal  or  pneumatic  leads  were  permitted  to  cross  the  bearing. 
This  latter  requirement  meant  that  dynamic  calibration  of 
the  system  would  be  complicated,  as  described  later,  but 
the  achievement  of  repeatable  mechanical  connections  be¬ 
tween  the  bearing  rotor  arid  stator  components  is  not  pos¬ 
sible  within  the  limits  acceptable  for  this  type  of  test¬ 
ing.  The  measurement  of  model  position  was  accomplished 
within  these  restraints  oy  the  design  and  fabrication  of 
a  unique  variable  reluctance  angle  transducer. 

Figures  4  and  5  show  the  bearing  and  model  assembly 
installed  in  the  wind  tunnel. 

The  model  is  fastened  to  the  air  bearing  rotor  by 
means  of  the  manifold  and  manifold  clamp  ring.  The  mani¬ 
fold  or  yoke  is  secured  to  the  rotor  and  distributes  the 
blowing  air  to  the  chambers.  The  angle  of  attack  trans¬ 
ducer  body  is  secured  to  the  air  bearing  housing  while 
its  armature  is  fastened  to  the  yoke. 

3.2  Mass  Transfer  System 

This  system  modulates  the  model  blowing  as  a  func¬ 
tion  of  the  angle  of  attack  and  includes  the  elements 
that  allow  calibration  so  that  instantaneous  values  of 
mass  flow  from  each  chamber  can  be  deduced.  Four  needle 
valves  modulate  the  gas  flowing  to  the  four  modal  com¬ 
partments.  Rate  of  mass  flow  is  inferred  from  four  small 
pressure  transducers  located  just  downstream  of  the  valves. 

The  four  model  compartments,  or  quadrants,  are  de¬ 
fined  by  the  intersection  of  the  horizontal  plane  through 
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the  longitudinal  axis  and  the  vertical  plane  that  contains 
the  axis  of  rotation.  Thus,  we  have  fore  and  aft  top  and 
fore  and  aft  bottom  compartments  which  are  shown  in  the 
sect i onal i zed  side  view  in  Figure  3A. 

The  path  of  the  biowing  air  from  the  needle  valve 
to  the  chamber  can  be  followed  by  reference  to  the  two 
assembly  views.  Upon  leaving  the  needle  valve  it  passes 
the  pressure  transducer  and  is  ducted  to  a  groove  on  the 
bearing  rotor  which  is  ducted  to  the  rotor  end.  The  yoke 
or  manifold  connects  the  four  rotor  outlets  to  their  re¬ 
spective  model  chambers. 

The  needle  valves  are  actuated  by  hydraulic  pressure 
to  close  and  pneumatic  to  open.  Lack  of  room  for  four 
return  lines  forced  the  use  of  a  common  pneumatic  return. 
Coupled  to  each  valve  stem  is  a  position  feedback  poten- 
t i ometer . 

The  pressure  transducers^  are  unique  because  of  their 
small  size  which  is  about  1  /4 11  O.D.  by  l/4"  long.  They 
aro  of  the  capacitance  type  and  ar e  used  in  a  special 
100  k.c,  carrier  system  which  us^s  triaxial  leads  to  cir¬ 
cumvent  changes  in  cable  capacitance.  The  relationships 
between  the  mass  flow  through  top  and  bottom  chambers  and 
the  angle  of  attack  are  given  in  Figure  6.  The  phase  lead 
shown  as  lambda  can  lag  as  well  as  lead.  The  fore  and  af 
sections  are  independent  of  each  other  and  can  be  set  at 
different  lambda  values.  The  valve  strokes  can  be  changed 
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independently  which  allows  adjustment  of  the  mass  flow 
ranges. 

The  functional  elements  of  this  system  for  one  valve 
are  shown  in  Figure  7.  The  needle  valves  are  positioned 
by  a  closed  loop  servo  system  which  is  programmed  by  the 
lead-lag  amplifier  which  has  the  angle  of  attack  system 
output  as  its  input.  The  hydraulic  servo  valves  were 
mounted  as  close  to  the  needle  valves  as  possible.  This 
meant  installing  them  at  the  base  of  the  model  support 
strut.  The  hydraulic  power  supply  was  outside  the  test 
section.  The  model  release  mechanism  allows  release  at 
+4  and  +8  degrees.  The  Hidyne  pressure  transducers  are 
mounted  on  the  inner  shell  of  the  model  with  one  side  of 
their  differential  input  in  the  blowing  compartments. 

Their  purpose  is  to  get  a  correlation  between  chamber 
pressure  and  the  pressure  back  at  the  Metrotech  pressure 
transducers  and  the  mass  flow  rate.  The  Hidyne  lead  wires 
are  removed  during  a  run.  For  dynamic  calibration  the 
signal  generator  simulates  the  model  motion.  The  static 
flow  calibration  system  contains  two  rotameter  type  flow 
meters,  dehumidifiers,  pressure  gage,  and  two  pressure 
regulators,  which  is  sufficient  equipment  to  obtain  the 
mass  flow  rate  at  steady  state  conditions.  Elements  of 
the  mass  transfer  control  and  calibration  system  are 
shown  on  Figures  8,  S,  10,  and  11. 

3.3  Models 


Two  types  of  models  were  used:  one  for  the  ablating 
series  and  another  for  the  blowing  series  of  runs.  Both 
types  used  the  same  inner  cone  shown  in  Figure  3A. 
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3.3.1  Ab  1  at i ns  Mode  1 

The  requirement  that  this  model  have  ablation  over 
the  whole  surface,  or  ablate  on  the  fore  or  aft  of  the 
CG  only,  leads  to  the  design  shown  in  Figure  12.  The 
coated  model,  when  initially  coated,  has  the  same  dimen¬ 
sions  as  the  blowing  model.  When  only  half  the  model  is 
coated,  the  other  half  is  replaced  with  a  similar  section 
of  greater  wall  thickness  to  form  a  smooth  juncture  be¬ 
tween  the  two  halves.  Two  solid  models  are  thus  possible: 
one  which  has  the  same  geometry  as  the  blowing  model  and 
another  which  !s  smaller.  Both  of  these  solid  models  were 
tested. 

The  model  geometry  is  as  follows:  (The  "bare"  fig¬ 
ures  are  the  dimensions  of  the  small  solid  model  and  the 
"coated"  figures  are  also  the  dimensions  of  the  large 
solid  model . ) 

A.  Cone  angle  -  10° 

B.  r/r.  =  .166 

n  b 

C.  Axis  of  rotation  =  56S»  of  length  from  the  nose 

D.  Open  base 

E.  Length  (bare)  =  11.03",  (coated)  =  11.65" 

F.  Base  diameter  (bare)  =  4.52",  (coated)  =  4.77" 

Ablation  models  are  shown  on  Figures  13,  14,  and  15. 

Two  different  materials  were  used  in  coating  the 
models  for  ablation  simulation  testing.  The  materials 
used,  parad  i  ch  lorobenzene,  C^H^C^,  and  ammonium  chloride, 
NH^Cl,  were  selected  because  they  each  sublime  at  the 
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temperature  levels  reached  and  because  they  differed  in 
mass  loss  rates.  Ammonium  chloride  was  applied  to  the 
nose  and  parad i ch 1 orobenzene  to  the  conical  section  of 
the  model. 

The  par  ad i ch  1  orobenzene  was  applied  to  the  model  in 
a  liquid  state  by  painting  it  onto  the  model  after  the 
model  surface  had  been  initially  prepared  by  coating  with 
Pliobond  cement.  When  a  sufficiently  thick  coating  of 
ablation  material  had  beer,  achieved,  the  model  was  ma¬ 
chined  to  the  desired  shape  with  the  use  of  metal  tem¬ 
plates.  The  resulting  surface  was  very  smooth  and  had 
a  waxy  texture. 

Prior  to  coating  the  nose  with  ammonium  chloride  it 
was  given  several  coats  of  flat  black  paint  and  sprinkled 
with  a  medium  coarse  grit.  This  surface  assured  a  good 
mechanical  bond  and  protected  the  mciel  from  the  corro¬ 
sive  effects  of  the  ablation  material.  The  finely  ground 
ammonium  chloride  used  was  mixed  with  water  to  form  a  thick 
paste  which  was  spread  on  the  model.  The  coated  model  was 
baked  until  the  moisture  was  removed  and  the  ammonium 
chloride  had  hardened.  The  coatings  were  machined  to 
shape  after  the  model  had  cooled.  Since  the  texture  of 
the  finished  surface  was  fine  grained  and  somewhat  crumbly 
the  model  was  handled  as  little  as  possible  after  coating. 

The  physical  properties  of  the  parad i ch 1 orobenzene 
and  ammonium  chloride  can  be  obtained  from  References  4, 

5,  and  6.  The  density  of  the  ablation  material  was  less 
than  theoretical  due  to  the  method  of  application.  Meas¬ 
urements  of  the  "as  applied"  material  density  have  been 
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3 

made  and  are:  1)  parad i ch 1 orobenzene  -  84.7  lbs/ft°  and 

3 

2)  ammonium  chloride  -  54.1  lbs/ft °. 

3.3.2  Bio.,  ng  Model 

The  design  objectives  of  this  model  were  that  the 
blowing  rate  decreases  for  increasing  downstream  stations 
and  that  the  air  volume  between  the  outer  porous  skin  and 
the  inner  shell  be  kept  to  a  minimum.  These  are  reflected 
in  Figure  3A  by  the  increasing  porous  skin  thickness  with 
increasing  station  and  tapering  of  the  blowing  compartments 
to  smaller  widths  away  from  the  gas  inlets.  By  keeping 
the  compartment  volumes  small,  their  response  is  increased. 

The  model  that  was  finally  tested  had  a  porous  skin 
aft  of  the  axis  of  rotation  and  a  solid  skin  forward  of 
it.  The  oorous  skin  was  Feltmetal^  and,  as  assembled, 
allowed  a  mass  flow  rate  of  approximately  .C004  lb/sec 
per  compartment  at  a  14  psia  internal  comparment  pres¬ 
sure  and  an  external  pressure  of  6  Torr.  The  blowing 
model  is  shown  on  Figure  5. 

3.4  Data  Acquisition  System 

The  requirements  of  this  system  are  that  it  record 
the  angle  of  attack  of  the  model  without  phase  lag  and 
render  its  amplitude  faithfully.  Also,  the  output  of  the 
pressure  transducers  from  which  the  mass  flow  rate  is  in¬ 
ferred  must  be  measured  with  high  fidelity.  The  model 
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oscillating  frequency  varied  from  2  *o  3.5  cps,  the  value 
for  any  wind  tunnel  run  being  a  function  of  the  stagna¬ 
tion  pressure  for  that  run. 

A  multi-channel,  light  beam  galvanometer,  oscillo¬ 
graph  was  used  to  record  all  the  data.  The  dynamic  data 
were  recorded  with  galvanometers  that  were  6456  of  criti¬ 
cally  damped  and  whose  natural  frequency  was  either  100 
cps  or  585  cps.  The  output/input  amplitude  ratio  for  both 
cases  is  one.  Tho  phase  lag  for  tha  100  cps  case  is  less 
than  1°  and  for  the  585  cps  case  is  approximately  0°.  The 
bandwidth  of  the  transducer  and  electronic  equipment  "feed¬ 
ing"  any  galvanometer  w?:  greater  than  the  galvanometer's. 
In  short,  the  data  acquisition  system  did  not  distort  the 
data. 
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4.0  OPERATING  PROCEDURES 


4.1  Ca  1  i  brat  i  o ns 


The  calibrations  necessary  are  the  bearing  damping, 
mode]  moment  of  inertia,  angle  of  attack  measuring  sys¬ 
tem,  static  and  dynamic  mass  transfer  measuring  systems, 
and  weight  loss/unit  time  of  ablative  material  in  the 
case  of  ablating  models. 

4.1.1  Bear i ng  Damp i ng 

This  is  determined  b>  clamping  a  gravity  pendulum 
to  the  model.  (See  Figure  16.)  This  combination  is 
allowed  to  oscillate  in  a  chamber  at  ambient  pressure 
and  a  fraction  thereof.  From  these  results  it  was  con¬ 
cluded  that  ambient  determination  can  be  extrapolated 
to  in  vacuo  conditions.  A  measure  of  damping,  the  log¬ 
arithmic  decrement  is  ottained  from  the  recording  of 
angle  of  attack. 

4.1.2  Moment  of  Inertia 


The  moment  of  inertia  of  the  oscillating  mass  was 
obtained  in  part  by  calculation  and  in  part  experimentally. 
Moment  cf  inertia  of  the  air  bearing  rotor,  oeing  of  a 
simple  shape  and  inaccessible,  was  calculated.  A  torsion 
pendulum  was  used  to  find  the  moment  of  inertia  of  the 
remaining  mass,  which  was  comprised  of  the  model,  mani¬ 
fold  ring  clamp  ana  manifold.  The  frequency  of  oscilla¬ 
tion  of  a  bar  of  known  moment  of  inertia  was  compared 
with  the  frequency  of  oscillation  of  •  unknown.  The 
apparatus  used  is  shown  in  Figure  17. 


14 


FluiDyne  Engineering  Corporal  i 


4.1.3  Angle  of  Attack 


This  measuring  and  recording  system  was  proven  with 
a  precision  inclinometer. 

4.1.4  Mass  Transfer 


The  calibration  of  the  mass  transfer  measurement 
system  was  carried  out  in  two  steps.  The  primary  stand¬ 
ard  employed  was  the  rotameter  which  has  a  dynamic  re¬ 
sponse  which  is  too  slow  to  be  used  to  calibrate  the  sys¬ 
tem  directly  under  dynamic  conditions.  The  procedure 
followed  was  to  use  the  rotameter  to  calibrate  dy, sarnie 
measuring  devices  (pressure  transducers)  under  steady 
flow  (static)  conditions,  after  which  a  dynamic  calibra¬ 
tion  was  made  using  on’y  the  pressure  transducers  to 
determine  the  lag  characteristics  of  the  transducers 
located  in  the  vicinity  of  the  flow  control  valves 
(which  are  monitored  during  a  test  run)  from  the  meas- 
i .  ement  of  the  pressure  in  the  chambers  by  the  chamber 
transducers  (which  are  not  monitored  during  a  test  run). 
These  ca  I  i  ations  were  made  with  the  test  cabin  evacu¬ 
ated  to  simulate  run  pressure  levels. 

4. 1.4. 1  Static  Calibration 

Tnis  was  done  by  metering  the  flow  to  one  compart¬ 
ment,  then  to  t^e  other  by  opening  and  closing  the  ap¬ 
propriate  needle  valves.  The  pressure  in  the  compart¬ 
ment  (Hidyne  transducer)  was  recorded  for  various  mass 
flow  rates.  Only  the  aft  top  and  bottom  compartments 
were  thus  calibrated  fc '  this  series  of  tests. 
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4. 1.4.2  Dynamic  Calibration 

Model  oscillation,  for  purposes  of  blowing  calibra¬ 
tion,  was  simulated  by  driving  the  valves  with  a  signal 
generator  at  the  frequency  expected  during  the  run.  Both 
the  compartment  pressure  and  the  pressure  at  the  needle 
valves  were  recorded.  This  proviued  a  correction  of  the 
pressure  at  the  needle  valve  with  the  flow  from  the  com¬ 
partment  through  the  previous  static  calibration. 

Compartment  cross-flow  was  checked  by  biasing  one 
valve  closed  and  observing  its  compartment  pressure 
fluctuations  es  the  other  valve  was  modulated.  Cross- 
flow  was  not  significant. 

4.2  Run  Procedures 


4.2.1  General 

The  operating  procedure  for  a  typical  run  was  as 
follows.  Prior  to  tunnel  pressurization,  the  model  was 
balanced  and  manually  set  to  one  of  the  four  detent 
positions  (+4°  and  +8°).  The  tunnel  was  then  closed 
and  pressurization  started.  After  reaching  steady  state 
P  ,  the  model  was  injected  and  the  detent  released.  Dur¬ 
ing  the  run,  the  time  history  of  the  model  oscillation 
was  recorded  on  an  oscillograph.  Prior  to  the  flow 
breakdown  the  model  was  retracted  from  the  stream. 

4.2.2  A  b  1  a  t  i  o  n  Runs 

For  these  runs  the  procedure  was  geared  to  obtain¬ 
ing  the  rate  of  sublimation  of  the  coating.  Time  inter¬ 
val  in  the  test  flow  was  ascertained  by  recording  model 
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support  position.  Upon  retraction  the  mode’  was  cooled 
by  a  jet  of  air  to  reduce  Siblimation  (see  Figure  15). 

For  a  typical  run  the  air  bearing  damping  was 
checked,  the  angle  of  attack  measuring  system  cali¬ 
brated,  the  model  weighed  and  mounted,  the  run  made, 
the  model  retracted  and  cooled  as  the  tunnel  wa»  shut 
down,  the  test  cabin  opened,  and  the  model  weighed. 

Moment  of  inertia  change  during  the  run  was  not  signi¬ 
ficant. 

4.2.3  Blowing  Runs 

Blowing  runs  differed  from  ablation  runs  only  with 
respect  to  the  differing  model  setup  and  handling  re¬ 
quirements.  Thus,  for  blowing  runs  the  principal  unique 
requirements  were  to  establish  the  blowing  characteristics 
desired  for  the  run,  and  to  calibrate  the  various  trans¬ 
ducers  involved  with  respect  to  the  data  acquisition  sys¬ 
tem  references  and  scales.  This  procedure,  while  simple 
in  principle,  requires  extreme  care  in  the  coordination 
of  the  several  subsystems  involved. 

Of  the  pressure  transducers  only  the  Hidyne  were 
accessible  for  calibration,  but  from  this  the  calibrations 
of  the  Metrotech  transducers  were  inferred.  Being  differ¬ 
ential-type  transducers,  pulling  a  vacuum  to  various  ievels 
on  the  side  not  in  the  model  blowing  compartment  gave  a 
calibration.  Then  by  exhausting  the  test  cabin  and  modu¬ 
lating  the  needle  valve  via  the  signal  generator,  a  Hidyne- 
Metrotech  correlation  was  achieved.  Also,  the  maximum  and 
minimum  mass  air  flows  from  each  compartment  were  set  and 
matched  by  adjusting  the  supply  pressure  and  valve  stroke 
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and  bias.  Then  the  test  cabin  was  opened  and  the  Hidyne 
wiring  removed . 

Again  the  air  bearing  damping  was  checked,  and  the 
angle  of  attack  measuring  system  calibrated  prior  to 
t ne  run,  as  in  the  abiation  runs. 
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5.0  DATA  REDUCTION  PROCEDURES  AND  RESULTS 


5.  1  Gfc.iera  1 


There  were  21  wind  tunnel  runs,  subdivided  as  follows: 

1.  four  runs  with  bare,  non-ablating  models  intended 
to  check  the  wind  tunnel  operating  characteris¬ 
tics  and  establish  the  aerodynamic  characteris¬ 
tics  of  the  mcdel  configuration  wi+hout  ablation 
or  blowing; 

2.  six  ablation  runs  involving  three  different  ab¬ 
lation  configurations;  and, 

3.  eleven  runs  with  the  blowing  model  of  which  one 
run  was  made  without  blowing,  for  reference  pur- 

•  poses,  and  the  other  ten  runs  involved  different 
phase  shift  angles  between  the  blowing  cycle  and 
the  model  oscillation  angle  and  different  blow¬ 
ing  intensities.  All  the  blowing  tests  were 
made  with  the  aft  blowing  configuration. 

Details  and  dimensions  of  the  various  models  are  given 
in  Section  3.3  and  in  Figures  3A  and  12.  Table  1  provides 
a  summary  of  the  runs  with  the  bare  and  ablating  models 
and  Table  2  summarizes  the  runs  with  the  biowing  model. 

The  runs  marked  with  an  asterisk  were  considered  to  pro¬ 
vide  valid  data  for  the  calculation  of  the  stability 
derivatives  and  these  runs  were  analyzed  in  detail. 

Data  from  the  other  runs  were  not  considered  valid  for 
the  reasons  summarized  in  the  Remarks  column  and  dis- 
c ussed  i n  Sect  ion  6. 
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information  pertinent  to  a  particular  run  and  neces¬ 
sary  for  the  data  reduction  is  extracted  from  the  follow¬ 
ing  sources. 

1.  The  run  log,  giving  the  tunnel  operating  con¬ 
ditions  PQ ,  T  ,  M,  and  the  model  geometric  and 
physical  characteristics. 

2.  The  oscillograph  traces  for  the  full  run  chart 
including  the  pre-calibrations  for  the  various 
transducers  and  the  full-run,  time-history  of 
the  model  angle  a,  test  cabin  static  p  essure 
p,  the  transducers  monitoring  P  ,  T  ,  or  in¬ 
dicating  model  injection  and  ejection,  and  for 
the  blowing  models,  the  pressures  in  both  top 
and  bottom  blowing  chambers  and  the  valve  posi¬ 
tion  transducer  outputs. 

3.  The  record  of  the  bearing  damping  check  carried 
out  prior  to  and  subsequent  to  the  run;  it  gives 
the  time  history  of  the  free  oscillation  of  the 
angle  a  for  the  pendulum  model  arrangement. 

A  photographic  record  complemented  the  data  for  each 
run.  The  ablation  runs  were  photographed  by  means  of  a 
schlieren  sequence  camera,  capable  of  shooting  up  to  10 
photographs  during  a  30-second  run.  Examples  of  the 
photographs  obtained  are  shown  on  Figure  18.  The  in¬ 
stants  of  exposure  are  recorded  by  identifiable  marks 
on  the  oscillograph  recording  chart.  A  continuous 
schlieren  photographic  record  was  made  for  all  the  blow¬ 
ing  :‘uns  using  a  16-millimeter  Mitchell  camera. 

5.2  Test  Section  Flow  Properties 


The  test  section  Mach  number  is  calculated  from  the 
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tunnel  calibration  curve  (Figure  2)  ?it  the  run  values 
of  and  T  .  Other  test  section  fiow  properties  are 
calculated  from  the  flow  tables  in  NACA  Report  1135 
with  the  Mach  number  as  the  independent  variable  and 
corrected  for  caloric  imperfections  by  means  of  the 
fo 1 lowi ng  equat ions. 


P  =  |P/P0KTP/P)p/p  (Po) 
q  =  (q/P0MTF/Plq/po(P0l 

t  =  <t/t0)(tp/p)t/t  (t0i 

0 

P  (P/P0)(TP/P)p/p 
0  '  o 

P  '^o  (T/T0)(TP/P)tAo 

v  =  49.  i  m yr 

where  (TP/P)  is  the  correction  factor  for  caloric  imper¬ 
fection  for  the  parameter  identified  by  the  related  sub¬ 
script.  The  factors  used  to  correct  for  caloric  imper¬ 
fections  are  obtained  from  plots  based  on  NACA  Report 
1135.  Also,  the  test  section  Reynolds  number  per  foot 
can  be  determined  by  a  su:*able  plot  based  on  the  same 
reference. 


5.  3  Stability  Derivatives 


The  static  and  dynamic  stability  derivatives,  M_ 

C  C  a 

and  M  +  M^  are  calculated  by  the  following  procedure. 


1.  First,  the  run  trace  is  carefully  examined  and 
any  aspect  of  irregular  behavior,  particularly 
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as  regards  the  cabin  static  pressure,  is  ex¬ 
amined  and  the  run  intervals  free  from  inegu- 
lsrities  are  selected. 

2.  The  test  section  flow  properties  are  obtained 
as  indicated  in  5.2. 

3.  The  galvanometer  deflections  for  amax  and 

a  .  are  read  out  from  the  oscillograph  trace 
min 

and  the  oscillation  frequencies  are  calculated. 

4.  From  the  a  calibration,  the  values  of  0  = 

1  /2 ( c*  w  -  a  .  )  are  calculated  and  the  time 
max  min 

history  curve,  0  versus  time,  is  plotted  on 
semi-  log  paper. 

5.  Choosing  suitable  time  intervals.  At,  on  the 
0-time  curve,  a  straight  line  tangent  to  the 
0  curve  is  drawn  for  each  time  increment  and, 
from  its  intercepts,  0j  and  0 at  the  ends 
of  the  interval,  the  quantity  1  n  ( 0  |/©2 )  J/At 

is  determined  as  an  approximation  to  the  deri¬ 
vative  -jjj  1  n  ( d  |/ ©2 )  • 

6.  The  total  moment  of  inertia  1^  ^  is  obtained  by 

summation  of  the  model  moment  of  inertia,  I  mod Jl 

and  the  bearing  moment  of  inertia,  I.  „. 

M*  Dea 

7.  The  total  damping  moment,  0^.^,  is  calculated 
from  the  equat  ion 


% 


tot 


=  -2(,tct) 


1  n  ( 0 ®  2 ) 


At 

M 


©bear *  ' s  °bta i ned 


8.  The  bearing  damping  moment, 

directly  from  the  bearing  damping  calibration 
and  the  value  of  t;.e  aerodynamic  damping  moment 
is  calculated  from  the  relation 


M0 

aero 


M; 


M; 


'tot 


bear 
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9.  The  static  stability  derivatives  are  calculated 
from  the  equat  ion 

.  2 

P  T  O  f  ^ 

- - — —  per  radian 

a  qSd 

10.  The  dynamic  stability  derivatives  are  calculattd 
from  the  relation 

Me 

CM  +  =  — ae-—  per  rad  i  an 

q  a  2 

qSd 

2V 

5.4  Mass  Transfer  Characteristics 


5.4.1  Ablation 

The  ablation  rate,  m,  in  IbsAec,  was  determined  by 
dividing  the  model  weight  loss  during  the  run  by  the  total 
run  duration.  In  Table  1,  the  ablation  intensity  is  also 
expressed  by  the  non-dimensional  quantity  ih/(pVS),  using 
the  free  stream  density  and  velocity  and  the  model  base 
area  for  reference  purposes.  It  should  be  pointed  out 
that  neither  m  nor  m/(pVS)  gives  a  true  indication  of 
the  local  or  even  the  average  rate  of  ablation  for  dif¬ 
ferent  ablation  configurations  of  the  same  model  since 
the  actual  ablation  area,  S^,  will  vary  for  the  same 
base  area.  More  realistic  parameters  for  comparing 
different  ablation  geometries  are  rti/ ( ^ }  and  m/lpVS^). 

The  parameters  m  and  m/(pVS)  are  used  only  for  the 
purpose  of  determining  the  corresponding  blowing  rates 
for  the  same  p°ometric  configuration. 


23 


FluiPyne  Engineering  Corporation 


5.4.2  J3j owi  ng 

One  of  the  principal  objectives  of  a  blowing  model 
is  to  simulate,  as  closely  as  possible,  t.he  mass  trans¬ 
fer  characteristics  of  the  ablating  body  over  the  full 
cycle  of  oscillation.  There  are  practical  limitations 
of  model  design  and  construction  which  make  exact  simu¬ 
lation  impossible  to  achieve.  In  the  present  model,  the 
continuous  transverse  variation  of  ablation  rate  from  a 
maximum  at  the  windward  side  to  a  minimum  at  the  leeward 
side  is  replaced  by  a  stepwise  change  between  a  top  and 
bottom  chamber.  The  downstream  variation  in  the  ablating 
mass  transfer  rate,  which  will  be  different  for  the  laminar 
and  turbulent  boundary  layer  cases,  is  approximated  ex¬ 
perimentally  by  increasing  the  skin  thickness  with  in¬ 
creasing  values  of  x.  Attempts  to  measure  the  distri¬ 
bution  of  mass  flow  through  the  porous  skin,  or  to  check 
the  actual  porosity  distribution  of  the  skin  by  means 
of  a  hot  wire  anemometer,  were  unsatisfactory  on  account 
of  the  very  small  flow  velocities.  Tests  of  the  blowing 
skin  under  water  suggested  the  presence  of  pinholes  and 
other  irregularities:  it  has  been  assumed  that  such  irregu¬ 
larities  were  randomly  distributed  and  do  not  unduly  viti¬ 
ate  the  experimental  observations. 

The  procedure  for  determining  the  timewise  distribu¬ 
tion  of  the  rate  of  air  mass  blowing  from  each  chamber  for 
a  particular  run  and  its  phase-shift  with  respect  to  the 
oscillation  angle,  a,  was  as  follows.  A  suitable  time 
interval  cf  2  cycles  of  oscillation  was  selected  from 
the  run  chart  and  divided  into  32  equal  intervals.  The 
galvanometer  traces  from  the  Metrotech  transducers  re¬ 
cording  the  blowing  chamber  pressures  were  read  off  and 
converted  to  mass  flow  in  lbs/sec  (by  use  of  the  combined 
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static-dynamic  calibrations  described  above)  for  each 
chamber.  The  mass  flow  was  plotted  along  with  the  a- 
trace  which,  for  all  practical  purposes,  was  a  sine 
curve.  The  resulting  mass  flow  distribution  was  not 
quite  sinusoidal  in  shape,  because  of  practical  experi¬ 
mental  limitations.  A  graphical  procedure  which,  in 
effect,  replaced  the  actual  waveform  by  an  appropriate 
sine  wave  with  the  same  amplitude  variation  was  adopted. 
This  enabled  an  average  phase-shift  angle,  assumed  con¬ 
stant  over  the  full  cycle,  to  be  determined.  In  general, 
the  mass  flow  rates  from  the  two  chambers  and  the  phase 
shift  angles  were  not  equal  in  spite  of  painstaking  set¬ 
tings  prior  to  the  wind  tunnel  run.  This  result  was 
due  to  the  fact  that  it  was  not  possible  to  set  up  the 
mass  flow  distributions  and  phase-shifts  with  the  actual 
model  in  oscillation  at  the  actual  test  frequency.  In¬ 
stead,  an  ct-simulator  was  used  for  the  initial  settings, 
with  the  model  itself  stationary.  The  addition  of  a 
torsion  pendulum  to  the  syster  would  eliminate  this  prob- 
1  em. 


It  is  worth  noting  that  a  theoretical  analysis  of 
a  blowing  model  may  be  employed  to  account  for  the 
effects  of  unequal  blowing  intensities  or  different 
phase  angles  on  the  static  and  dynamic  stability  deri- 
vat i ves. 
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5. C  -NALYSIS  OF  RESULTS  AND  CONCLUSIONS 


6. 1  General 


As  indicated  in  Section  5,  tl.e  data  from  some  runs 
had  to  be  discarded.  The  main  cause  for  invalidating  a 
test  run  was  a  persistent  cabin  static  pressure  oscilla¬ 
tion  of  the  same  frequency  as  that  of  the  model  oscilla¬ 
tion;  when  this  pressure  fluctuation  was  appreciable,  it 
appeared  to  drastically  affect  the  ralculated  stability 
der i vat i ves . 

Another  type  of  irregularity,  observed  with  both  the 
ablation  and  blowing  models,  was  that  of  random,  rela¬ 
tively  large-amplitude,  cabin  static  pressure  pulses  which 
had  a  strong,  though  localised,  effect  on  the  model,  readily 
evident  in  the  abrupt  change  in  the  a-amplitude.  Initially, 
it  ,uas  feared  that  these  static  pressure  transients  viti¬ 
ated  the  data  for  the  entire  run.  However,  careful  exam¬ 
ination  of  the  records  of  a  number  of  runs  showed  satis¬ 
factory  consistency  of  the  results  calculated  from  the 
data  for  the  intervening  periods  of  steady  cabin  static 
pressure,  particularly  when  these  quiescent  periods  were 
of  relatively  long  duration,  of  the  order  of  5  seconds 
or  longer. 

The  plot  of  the  time  history  for  the  run  (0  versus  t) 
on  semi  log  paper  was  found  to  be  a  very  convenient  refer¬ 
ence  for  an  initial  analysis  of  the  run.  Tha  intervals 
o'  time  during  which  the  data  are  not  influenced  by 
static  pressure  transients  (arid  thus  are  most  suitable 
for  reduction/  are  clearly  shown. 
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6 .  2  Ablation  Tes t_s 


6.2.1  Air  Bearing  re.  fprmance 

The  air  bearing  performed  quite  well  for  these  runs 

the  bearing  damping  being  small  and  generally  repeatable 

A  plot  of  bearing  damping  is  shown  in  Figure  19.  The 

general  level  of  damping  appeareH  constant  with  ampli- 

tude  and  a  value  of  .05  x  10  f t- 1  b/r ad-sec  was  used 

for  the  data  reduction  of  these  runs.  For  comparison, 

the  average  damping  calculatei  for  the  earlier,  strut- 

-4  / 

mounted  air  bearing  was  .22  x  10  f t- 1  b/ rad-sec . 

It  should  be  pointed  out  that  the  gravity  pendulum 
used  for  the  bearing  damping  checks  allowed  the  system 
to  oscillate  at  approximately  1  cps,  whereas  the  fre¬ 
quency  observed  during  the  actual  runs  ranged  between 
2  and  3.5  cps . 

6.2.2  Spec  i  f i ;  Remark s  and  Resu  Its  of  t he 
Ablation  Test  Runs 


Taole  1  presents  a  convenient  summary  of  the  abla¬ 
tion  test  results,  indicating  the  successful  runs,  and 
Figure  20  shows  the  time  history,  0  versus  t,  foi  these 
runs.  The  following  remarks  are  pertinent  to  specific 
runs  of  Table  1. 

Run  845,  with  the  small  bare  model,  involved  only 
slight  cabin  pressure  oscillations  and  no  transient 
static  pressure  pulses.  The  model  was  released  at  8 
degrees,  ^ s c i 1 1  a f e d  at  a  frequency  of  about  2.5  cps, 
and  converged  as  expected.  Table  3A  lists  some  of  the 
relevant  characteristics  for  this  run  together  with  the 
calculated  static  and  dynamic  derivatives. 
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Run  846  involved  the  large  model  in  the  aft-ablation 
configuration.  The  model  was  released  at  4  degrees  and 
its  oscillations  diverged.  The  parad i ch  lorobenzene  cub- 
limateH  at  the  rate  of  .006  lb/sec  and  slight  oscillations 
were  observed  in  the  traces  of  the  cabin  static  pressure. 
Details  of  this  run  and  the  calculated  derivatives  are 
listed  in  Table  3B. 

Rut  847  was  essentially  a  repetition  of  Run  845. 
However,  it  displayed  a  number  •  ,  random,  static  pres¬ 
sure  pulses  during  ine  initial  seconds  of  the  run,  after 
which  period  the  cabin  static  pressure  remained  fairly 
steady  and  the  subsequent  data  reduction  yielded  deri¬ 
vatives  compatible  with  those  of  Run  845,  thus  support¬ 
ing  the  assumption  that  valid  results  may  be  calculated 
from  data  away  from  the  local  static  pressure  disturb- 
nce.  Table  3C  presents  the  results  obtained  from  this 
run . 


Run  848  with  the  large,  f u 1 1- 1 engt h-ab  lat i on  model, 
and  Runs  849  a*»d  850  with  the  same  model  in  the  forward- 
ablation  configuration,  all  displayed  relatively  large, 
persistent  oscillations  in  the  cabin  static  pressure,  at 
the  same  frequency  as  the  model  oscillation,  resulting  in 
divergent  oscillations  for  all  three  runs.  It  was  sus¬ 
pected  that  the  model  was  marginally  large  for  the  tunnel 

operat:on  conditions  of  P  =  500  psia  and  T  =  3200°R. 

o  o 

Run  852,  conducted  with  a  bai  e  model  the  same  size 
as  the  ablating  models,  at  PQ  =  500  psia  and  TQ  =-  3100°R, 
confirmed  the  above  suspicions.  Appreciable  oscillations 
in  the  cabin  static  pressure  were  observed  and,  contrary 
to  the  expected  results,  as  demonstrated  by  Runs  845  and 
847,  the  model  oscillations  diverged. 
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The  retraining  runs  in  this  series,  nambely  851,  853, 
and  854,  were  made  at  approximately  PQ  =  1000  psia  and 
T0  =  2500°R  and  they  all  exhibited  a  steady  cabin  static 
pressure  and  produced  valid  run  data.  Test  Runs  851  and 
853  with  the  large,  forward  ablation  configuration  and 
release  angles  of  4  and  8  degrees  respectively,  dis¬ 
played  the  same  degree  of  convergence  and  approximately 
the  same  values  for  the  derivatives  as  shown  in  Tables 
3D  and  3E.  The  results  for  the  final  test  run  in  this 
series.  Run  854  with  the  small, Dare  model  released  at 
8  degrees,  are  given  in  Table  3F.  The  calculated  values 
for  the  static  and  dynamic  derivatives  are  compatible 
with  the  values  obtained  from  Runs  845  and  847  at  differ¬ 
ent  wind  tunnel  operating  conditions. 

Figures  21  and  22  illustrate,  for  purposes  of  com¬ 
parison,  the  calculated  static  and  dynamic  derivatives 
*rom  all  the  valid  runs  in  this  series. 

Figure  21  shows  that,  for  the  CG  location  of  the 

present  models,  the  variations  in  the  static  derivative, 
r 

Ma,  with  forward  aft  ablation  are  relatively  small, 

of  the  order  of  10  percent  of  the  original  value  of 
r 

Ma.  This  variation  was  stabilizing  in  the  forward 
ablation  configuration  anc  destabilizing  in  the  aft 
ablation  configuration. 

Figure  22  illustrates  the  variations  in  the  dynamic 
C  C 

derivative,  M  +  M*,  produced  as  a  result  of  forward 

Q 

or  aft  ablation.  The  variation  in  amplitude  is  roughly 
one  order  of  magnitude  larger  than  the  original  value  and 
is  stabilizing  for  the  forward-ablation  and  destab' 1 i zi ng 
for  the  aft-ablation  configurations. 
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Figure  23  i  1  lustrat «.  j  a  typical  calibration  curve 
of  the  a-transducer  takpn  from  one  of  the  ablation  tests. 

6. 3  B lowi ng  Tests 


6.3.1  Air  Bearing  Performance 


The  air  bearing  performance,  as  indicated  by  the  meas- 

M' 

ured  bearing  damping  moment,  0^  ,  's  'l^^rated  in 

Figure  24  for  the  valid  blowing  test  runs.  Straight  line 
approximations  to  the  test  values  are  also  shown.  For 
data  reduction  purposes,  the  value  of  the  'earing  damp¬ 
ing  for  each  test,  at  the  relevant  angle  of  oscillation, 
was  read  off  the  appropriate  straight  line.  It  is  appar¬ 
ent  that,  even  for  the  best  damping  case,  associated 
with  test  Run  881,  when  the  bearing  damping  moment  was 
the  least  and  very  nearlv  constant  for  the  whole  range 
of  oscillation  angles,  its  numerical  value  of  about 
.4  x  10  f t- 1 b/rad-sec  was  considerably  larger  than 

_4 

that  of  the  ablation  runs  which  averaged  only  .05  x  1C 
f  t -  1  b/ rad-sec. 


All  the  bearing  damping  tests  were  cart  i ed  out  with 
a  pendulum  arrangement  which  oscillated  the  model  at  a 
frequency  close  to  1  cps;  whereas  the  frequency  of  os¬ 
cillation  for  the  blowing  tests  centred  around  2.5  cps. 

With  the  larger  bearing  damping  moments  associated 
with  the  blowing  model,  there  was  no  appreciable  change 
in  the  measured  bearing  damping  when  the  damping  tests 
were  carried  out  in  the  atmosphere  or  in  an  evacuated 
chamber. 
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The  design  demands  on  the  air  bearing  for  the  multi- 
chamber,  blowing  model  are  far  more  stringent  than  those 
for  the  ablation  model.  In  the  latter,  the  air  bearing 
essentially  has  only  to  provide  a  free  support  for  the 
model  oscillation;  whereas  in  the  former  it  has  also  to 
channel  and  control  the  blowing  air  to  the  respective 
chambers.  Some  difficulties  with  bearing  heating  and 
overloading  were  experienced  when  attempting  to  run  the 
blowing  tests  at  the  wind  tunnel  conditions  which  had 
given  consistently  steady  cabin  static  pressures  for  the 
ablation  runs,  namely  PQ  =  1000  psia  and  Tq  =  2500°R. 

6.3.2  Specific  Remarks  and  Results  of  the 
Blowing  Test  Runs 


Table  2  presents  a  summary  of  the  11  blowing  test 
runs  and  indicates  with  an  asterisk,  the  7  runs  which 
appeared  to  provide  valid  data  suitable  for  detailed 
reduction  and  analysis. 

The  wind  tunnel  conditions  chosen  for  Test  Runs  874, 
875,  and  876  were  those  found  to  be  the  best  for  the 
ablation  tests.  The  release  angle  of  4  degrees  for  these 
runs  was  chosen  in  anticipation  of  diverging  oscillations 
for  an  aft  blowing  configuration  with  phase  iag.  How¬ 
ever  these  three  runs  showed  very  marked  • rr egu  lar i t i es 
and  had  to  be  discarded  Run  874  exhibited  the  oscillat¬ 
ing  cabin  static  pressure  described  previously  with  rap¬ 
idly  decaying  oscillations  at  a  frequency  of  about  3.5 
cps . 


Run  875  was  plagued  witn  a  delayed  operation  of  the 
valve  in  the  top  chamber,  an  increasing  model  ce/ity  pres¬ 
sure,  and  very  rapidly  decaying  model  oscillations. 
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A  large,  initial  static  pressure  pulse  occurred 
near  the  model  release  for  Run  876  and  it  oscillated  for 
barely  two  cycles.  After  each  of  the  above  three  runs, 
the  model  became  very  hot  and  noticeably  sluggish.  Damp¬ 
ing  checks,  performed  immediately  after  the  runs,  indi¬ 
cated  a  considerable  deterioration  in  the  bearing  per¬ 
formance,  but  the  bearing  damping  characteristics  re¬ 
turned  to  their  earlier  level  after  the  model  had  re¬ 
turned  to  room  temperature  conditions. 

All  subsequent  runs,  namely  877  to  884,  were  made 

at  around  500  psia  for  P  ,  and  2400°R  for  T  ,  and  th:s 

0  0 

alleviated  the  force  and  thermal  loading  on  the  bearing 
and  model  which  remained  relatively  free  and  only  moder¬ 
ately  hot  after  these  runs.  The  general  level  of  the 
bearing  damping  crept  up,  however,  with  successive  runs 
as  s hown  in  Pi gur e  24. 

In  view  of  th  •  fact  (already  noted  for  the  ablation 
runs)  that  the  wind  tunnel  operating  conditions,  in  terms 
of  diffuser  performance,  were  somewhat  marginal  for  the 
present  size  of  the  model,  random  static  pressure  pulses 
invariably  occurred,  sometimes  in  large  number  as  indi¬ 
cated  in  Table  2  and  Figure  25.  This  tended  to  restrict 
the  valid  data  from  some  runs  to  shorter  intervals  and 
may  be  the  cause  of  greater  scatter  in  the  calculated 
der i vat i ves. 

Table  2  lists  the  maximum  and  minimum  blowing  rates 
per  chamber  and  the  phase-shift  angles  X^  and  X2  for  the 
upper  and  lower  chambers  respectively,  calculated  by  the 
method  uilined  in  Section  5.4.2. 
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Improper  cycling  of  the  upper  chamber  valve  some¬ 
times  occurred  and  it  remained  nearly  fully  open  (Run 
877)  or  fully  closed  (Run  882)  and  thus  did  not  contrib¬ 
ute  to  the  derivative.  The  phase-shift  angles  for  the 
upper  and  lower  chambers  differed  markedly  in  the  early 
runs.  In  the  last  two  runs,  however,  as  a  result  of 
working  experience  in  setup  and  calibration,  the  phase- 
shift  angles  for  the  two  chambers  were  remarkably  close. 

Tables  4A  to  4G  list  the  essential  results  of  all 
the  valid  b lowi ng  r uns . 

The  calculated  values  of  the  static  stability  deriv- 
r 

ative,  Ma,  for  the  bluwing  runs  are  collected  together 

and  illustrated  in  Figure  26  which  also  shows  the  range 
r 

of  values  of  M  for  the  non-ablating  models  to  serve  as 

u 

reference.  It  will  be  noted  that  the  range  of  variation 
r 

in  Mj  due  to  blowing  at  the  different  rates  and  phase- 

u 

shift  angles  used  is  relatively  small,  of  the  order  of 
10  percent.  For  the  predominant  cases  of  phase-lag,  the 
effect  of  aft  blowing  is  seen  to  be  statically  destabil¬ 
izing. 

All  the  calculated  values  of  the  dynamic  stability 
C  C 

derivative,  M  +  M-  for  the  blowing  runs  are  collected 

q  a  C  C 

together  in  Figure  27.  The  range  of  variation  of  M  +  M^ 

H 

for  the  bare,  non-ablating  model  is  also  illustrated  for 
reference  purposes.  In  spite  of  the  large  scatter,  evidenced 
by  the  results  of  some  of  the  runs,  and  the  very  small  rates 
of  blowing  in  those  runs.,  the  following  trends  clearly  emerge 

1.  Run  878,  involving  negligible  phase-shift,  ex¬ 
hibits  no  appreciable  change  in  the  dynamic  deriv¬ 
ative. 
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2.  Runs  879  and  882  produced  the  maximum  positive 
dynamic  derivatives;  in  other  words,  the  effect 
of  blowing  was  strongly  destabilizing.  These 
runs  involved  the  largest  phase-lag  angles, 
although  in  Run  882,  only  the  lower  chamber 
cycled  satisfactorily 

3.  Run  883,  the  only  run  with  large  phase-lead 
angles,  was  the  only  run  where  the  blowing  in¬ 
creased  the  dynamic  damping  derivative. 

6. 4  Concluding  Remarks 


In  any  test  program  there  are  a  great  many  factors 
which  influence  the  precision  of  the  test  results.  In 
a  program  in  which  the  major  effort  has  been  the  develop¬ 
ment  of  the  test  hardware  and  a  test  technique,  ar  error 
analysis  becomes  impractical;  such  is  the  case  here.  The 
following  commentary  is  therefore  intended  to  give  a  basis 
for  evaluating  the  potential  capability  of  the  test  hard¬ 
ware  and  technique,  rather  than  make  any  attempt  to  pro¬ 
vide  a  classical  error  analysis  of  the  component  and  sys- 
ror  characteristics. 

The  test  facility  itself,  the  Hypersonic  Wind  Tunnel, 
provides  a  reliable  test  stream  in  which  to  evaluate  hyper¬ 
sonic  flow  characteristics.  The  run  time  of  the  order  of 
1  minute  permits  the  use  of  instrumentation  which  has  been 
widely  accepted  and  which  car  be  readily  calibrated.  One 
of  the  principal  factors  influencing  the  precision  of  co¬ 
efficients  determined  by  hypersonic  testing  is  the  prob¬ 
lem  of  determination  of  the  true  flow  Mach  number,  and 
thus  dynamic  pressure,  to  be  used  in  reducing  the  data 
to  coefficient  form.  In  tests  where  large  transverse  and 
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axial  gradients  occur,  the  absolute  values  of  coefficients 
can  become  meaningless.  The  use  of  a  contoured  nozzle 
in  the  FluiDyne  wind  tunnel  minimizes  the  latter  considera¬ 
tions  such  that  the  total  probable  variation  in  Mach  number 
over  the  test  core  and  the  error  in  its  average  level  can 
be  expected  to  be  within  1%.  That  this  is  the  case  is 
born  out  by  the  consistent  results  obtained  from  tests 
in  this  facility. 

The  special  test  apparatus  has  the  same  basic  char¬ 
acteristic  as  the  facility;  the  instrumentation  used  is 
of  a  type  which  has  been  generally  widely  accepted  and, 
in  this  case  more  importantly,  can  be  readily  calibrated. 

In  order  to  obtain  a  meaningful  correlation  between  mass 
flow  addition  characteristics  and  resultant  vehicle  sta¬ 
bility  characteristics,  it  is  necessary  to  know  what  the 
former  are  with  an  acceptable  degree  of  precision.  This 
is  the  basis  for  using  a  blowing  system  as  opposed  to 
employing  natural  ablation  or  sublimation  of  a  coating, 
i.e.,  to  achieve  control  over  the  experiment.  The  two 
principal  parameters  of  interest  are  the  blowing  rate 
and  its  phase  relationship  to  the  body  oscillation. 

Neither  of  these  parameters  can  be  precisely  determined 
for  the  ablation  or  sublimation  test.  In  the  case  of 
blowing  tests,  it  is  conceivable  that  these  parameters 
might  be  measured  directly;  however,  practical  con¬ 
siderations  generally  require  that  the  measurements  be 
made  indirectly  during  a  run  by  use  of  pre-run  calibra¬ 
tions  and  direct  measurement  of  related  parameters  dur¬ 
ing  a  run.  Mass  flow  calibrations,  in  effect,  calibra¬ 
tion  of  the  flow  through  the  porous  skin,  are  made 
statically,  i.e.,  with  constant  pressure  differential 
across  the  wall.  If  pressure  is  measured  upstream  of 
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the  chamber  during  a  run,  the  chamber  pressure  must  be 
inferred  from  calculations  or  a  dynamic  calibration  made 
before  the  run.  If  the  frequency  at  which  the  system  is 
operating  is  high  enough,  both  the  amplitude  decrement 
and  phase  shift  between  the  pressure  measuring  point  and 
the  chamber  can  be  substantial.  Thus,  although  the  flow 
path  between  the  measuring  point  and  the  mass  injection 
point  may  be  short  (such  that  the  time  required  for  a 
pressure  pulse  to  traverse  it  is  negligible),  the  impor¬ 
tant  consideration  is  that  of  the  flow  process  which  occurs 
between  these  two  points.  For  the  system  described  here 
which  operates  at  a  frequency  of  several  cycles  per  sec¬ 
ond,  dynamic  calibrations  have  demonstrated  that  the 
effect  on  amplitude  and  phase  shift  is  negligible.  (In 
the  case  where  the  system  frequency  is  an  order  of  magni¬ 
tude  larger,  both  of  these  factors  can  be  expected  to  be 
significant,  and  extremely  difficult  if  not,  practically 
speaking,  impossible  to  calibrate  out.) 

In  the  tests  described  herein  it  was  found  that  even 
with  the  low  frequency  system  and  the  provision  of  means 
for  dynamic  calibration,  difficulties  in  pre-run  setup 
and  post-run  data  analysis  were  encountered.  These  were 
largely  due  to  the  fact  that  the  dynamic  calibration  was 
achieved  by  use  of  a  stationary  model  and  a  function  gen¬ 
erator  to  provide  the  input  to  the  flow  control  system, 
or  an  oscillating  model  for  which  a  gravity  pendulum  pro¬ 
vided  the  driving  force.  In  the  latter  case,  the  fre¬ 
quency  of  oscillation  was  only  a  fractior  of  that  during 
wind  tunnel  testing;  thus,  both  approaches  only  achieved 
partial  simulation  of  the  complete  dynamic  system.  Any 
further  testing  should  incorporate  the  use  of  a  torsion 
pendulum  as  a  device  to  oscillate  the  model  at  test  fre¬ 
quencies  with  all  systems  operational  to  provide  final 
dynamic  calibrations  and  pre-run  setup. 
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The  single  technical  problem  which  gave  rise  to  the 
largest  number  of  practical  difficulties  in  the  experi¬ 
mental  program  was  locating  a  suitable  material  for  the 
porous  surfaces  of  the  blowing  model.  A  suitable  mate¬ 
rial  for  this  purpose  should  have  the  following  proper¬ 
ties: 

1.  Fine-grained  porosity  so  that  when  gas  is  blown 
through  it,  the  length  scale  of  nonuniformities 
in  local  mass  flow  rate  will  be  small  compared 
to  the  boundary  layer  thickness. 

2.  A  porosity  high  enough  that  the  permeability  of 
practical  thicknesses  is  high  enough  to  pass 
the  desired  flow  rates  at  chamber  pressure  low 
enough  that  the  metering  valves  remain  choked. 

3.  A  porosity  low  enough  that  the  blowing  chamber 
pressure  be  high  enough  so  that  the  volumetric 
flow  rate  between  the  control  valves  and  the 
chamber  is  low  enough  to  preclude  choking  be¬ 
tween  these  points. 

4.  Mach i nabi 1 i t y  or  formability  so  that  an  accurate 
cone  frustum  may  be  produced  without  closing  the 
pores  at  the  surface. 

5.  Sufficient,  mechanical  strength  to  withstand  the 
stresses  developed  by  the  internal  (chamber) 
pressure. 

6.  Uniformity  such  that  variations  in  porosity 
(averaged  over  areas  small  with  respect  to  the 
total  porous  area)  will  be  less  than  20%. 

During  the  course  of  a  year,  several  different  blow¬ 
ing  surfaces  were  fabricated  of  s i ntered-meta  1  products 
which  were  expected  to  have  these  characteristics;  none 
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of  the  resulting  surfaces  were  satisfactory  when  finally 
tested.  The  material  v/h'ch  gave  the  best  results  and 
which  was  used  in  the  test  program  reported  here  turned 
out  to  be  very  difficult  to  braze.  An  intermittent  crack 
developed  along  the  brazed  longitudinal  seams,  whi^h  were 
located  at  +90°  and  -90°  from  the  windward  surface  at 
positive  angle  of  attack  and  undoubtedly  a  significant 
amount  of  the  regulated  blown  air  supply  leaked  out  of 
these  cracks  instead  of  bleeding  through  the  surface. 

The  interaction  of  this  leaked  airflow  with  the  external 
hypersonic  flow  could  not  cause  much  change  in  pitching 
moment,  since  the  pressure  loads  thus  generated  would  be 
at  right  angles  to  the  pitching  plane.  It  is  thus  prob¬ 
able  that  the  actual  mass  transfer  which  was  effective: 
in  changing  the  pitching  moment  and  the  p:tch  dampinr 
rate  was  much  smaller  than  the  values  calculated  and 
tabulated  in  Table  2.  Also,  as  noted  previously,  there 
was  evidence  that  the  porosity  distribution  was  quite 
nonuniform  so  that  a  significant  amount  of  the  regulated 
blown  air  supply  may  have  been  blown  through  the  surface 
in  small  localized  regions.  This  would  also  tend  to  re¬ 
duce  the  effectiveness  of  the  blown  air  in  changing  the 
aerodynamic  pitching  moment. 

After  these  tests  were  completed,  we  located  a  third 
supplier  whose  sintered  metal  product  appears  to  satisfy 
the  practical  requiiements  listed  above.  We  plan  to  make 
bench  tests  of  a  model  fab.  cated  of  this  material  to 
determine  whether  this  is  the  case. 

The  limited  experimental  data  obtained  thus  far  and 
reported  herein  qualitatively  confirm  the  results  ob¬ 
tained  from  the  theoretical  analysis  of  the  phenomenon 
(Reference  3).  As  mentioned  previously,  that  portion 
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of  Ron  879  which  was  free  of  test  cabin  pressure  fluctu¬ 
ations  clearly  demonstrated  that  the  model  was  dynamically 
unstable  with  aft  blowing  with  phase  lag,  as  predicted 
by  the  analysis.  Similarly',  Run  883  showed  conclusively 
that  the  pitch  damping  rate  was  greatly  increased  for 

aft  blowing  with  phase  lead.  The  changes  in  t  e  dynamic 

C  C 

stability  derivatives,  ,  for  these  runs  was 

only  about  one-tenth  of  the  change  predicted  by  the  anal¬ 
ysis,  however.  A  large  part  of  this  quantitative  differ¬ 
ence  was  undoubtedly  the  result  of  the  ieaks  and  non- 
uniformities  in  the  porous  skin  described  above. 


( 
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